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Abstract 


The theoretical understanding of emergent phenomena in quantum materials is one of the greatest 
challenges in condensed matter physics. In contrast to simple materials such as noble metals and 
semiconductors, macroscopic properties of quantum materials cannot be predicted by the properties 
of individual electrons. One of the examples of scientific importance is strongly correlated electron 
system. Neither localized nor itinerant behaviors of electrons in partially-hlled Sd, 4/, and 5/ 
orbitals give rise to rich physics such as Mott insulators, high-temperature superconductors, and 
superior thermoelectricity, but hinder quantitative understanding of low-lying excitation spectrum. 
Here, we present a new first-principles approach to strongly correlated solids. It is based on a 
combination of the quasiparticle self-consistent GW approximation and the Dynamical Mean Field 
Theory (DMFT). The sole input in this method is the projector to the set of correlated orbitals for 
which all local Feynman graphs are being evaluated. For that purpose, we choose very localized 
quasiatomic orbitals spanning large energy window, which contains most strongly-hybridized bands 
as well as upper and lower Hubbard bands. The self-consistency is carried out on the Matsubara 
axis. This method enables the first-principles study of Mott insulators in both their paramagnetic 
(PM) and antiferromagnetic (AFM) phases. We illustrate the method on the archetypical charge 
transfer correlated insulators La 2 Cu 04 and NiO, and obtain spectral properties and magnetic 
moments in good agreement with experiments. 
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Introduction. The first-principles description of strongly-correlated materials is currently 
regarded as one of the greatest challenges in condensed matter physics. The interplay be¬ 
tween correlated electrons in open d- or /- shell and itinerant band states gives rise to rich 
physics that makes these materials attractive for a wide range of applications such as oxide 
electronics, high-temperature superconductors and spintronic devices. Various theoretical 


approaches are currently being pursued 


l|. One of the most successful approaches is the 


^namical mean held theory (DMFT) 2(]. In combination with density functional theory 
3|, 1^, it has described many features of strongly-correlated materials successfully and high¬ 
lighted the surprising accuracy of treating correlations local to a small subset of orbitals 
exactly, while treating the reminder of the problem in a static mean held manner. ^ . 

The numerous successes of DMFT in diherent classes of correlated materials revived 
the interest in the long-sought goal of achieving a diagrammatically controlled approach to 
the quantum many-body problem of solids, starting from the Green’s function G and the 


screened Coulomb interactions W j^, 7|. The lowest order diagram in perturbation theory 
in this functional gives rise to the GW approximation js] while the local approximation 
applied to the most correlated orbitals gives rise to an extended DMFT approach to the 
electronic structure problem l7|. The addition of the GW and DMFT graphs was proposed 

n n 

and implemented in model Hamiltonian studies [9[ and in realistic electronic structure [10, 
111 ]. There is now intense activity in this area with many recent publications 121-115| triggered 
by advances in the quality of the impurity solvers [l^, 17|, insights into the analytic form of 
the high-frequency behavior of the self-energy and improved electronic structure codes. 

Several conceptual issues remain to be clarihed before the long sought goal of a robust 
electronic structure method for solids is attained. The hrst issue is the choice of local 
orbitals on which to perform the DMFT method (summation of all local Feynman graphs). 
The second issue is the level of self-consistency that should be used in the calculation of 
various parts of the diagrams included in the treatment (free or bare Green’s function Go vs 
self-consistent interacting Green’s functions G). These central issues are addressed in this 
letter. 

The self-consistency issue appears already at the lowest order, namely, the GW level, and 
it has been debated over time. The corresponding issue in GW-I-DMFT is expected to be 


3 , 


20 


at least as important, but has not been explored, except for model Hamiltonians 
At the GW level, it is now well established that Hedin’s fully self-consistent formulation [Sj], 
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while producing good total energies in solids [2l|, atoms and molecules [22|, [2^, does not 


even 3D electron gas and aluminum in 
24| . Instead, using a free (quasiparticle) 


produce a good approximation to the spectra o: 
comparison to non-self-consistent GW results 
Green’s function in the evaluation of the polarization graph of the GW method gives much 
better results for spectral functions. This is the basis of the one-shot quasiparticle (QP) 
GW, starting from LDA 2^ or from others 2^. Unfortunately, the answer depends on 


the starting point. A solution for this problem is to impose a self-consistency equation 
to determine Gq. This method, called the quasiparticle self-consistent GW (QSGW) j^. 


is very successful reproducing the spectra of many systems 


DMFT is an important open challenge 


27l |. How to combine it with 


28|,l29|. 


Previous GW-I-DMFT studies typically used a Gq which depends on the LDA starting 


point, and projectors spanning a relatively small energy window 


12 


15|. In this work, we 


propose a different approach to the level of self-consistency and the choice of the DMFT 
orbital. We do a self-consistent QSGW calculation and then calculate local self-energy using 
DMFT with static Ud and Jh without feedback to non-local self-energy within GW. For the 
DMFT step, we choose a very localized orbital spanning large energy window which contains 
most strongly-hybridized bands as well as upper and lower Hubbard bands. 

In the LDA-I-DMFT context, the choice of very localized orbitals has provided a great 
deal of universality since the interactions do not vary much among compounds of the same 
family. This has been demonstrated in the studies of iron pnictides 30|] and transition metal 
oxides Ml- This choice results in a second advantage as we will show below, namely the 


frequency dependence of the interaction matrix can be safely ignored. Having chosen the 
correlated orbitals, all the other parameters are self-consistently determined. This is the 
first ab initio quasiparticle self-consistent GW-I-DMFT implementation and the first study 
on a paramagnetic Mott insulator within the GW-I-DMFT method. 

Results. Fig. 2(a) shows the frequency dependence of real and imaginary parts of Ud of 
La 2 Cu 04 shown in Fig. 1. It is calculated on an imaginary frequency axis and analytically 
continued by a maximum entropy method 32|]. We also plot the fully screened Coulomb 
interaction Wd for comparison. Static Ud is 12.0 eV and Ud remains almost constant up to 
10 eV. In contrast, in Wd, there are several peaks due to low-energy collective excitations 
below 10 eV. At very high energy, Ud approaches the bare coulomb interaction of 28 eV. 
Static value of Upd is 2.0 eV, much smaller than Ud, hence we don’t discuss its treatment 
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further 


33| . Calculated Jh is 1.4 eV and has negligible frequency dependence. By contrast, 


conventional constrained-RPA, in which 10 bands of mostly Cu-3d character are excluded 
from screening, results in static Ud = 7.6 eV, which is too small to open the Mott gap, and 
which is also inconsistent with photoemission experiments on CuO charge transfer insulators 


34|. 


We introduced a complementary method to compute the static Ud- The key idea is to 
hrst calculate the excitation spectra of La 2 Cu 04 within MQSGW+DMFT using local GW 
(with a static Ud) as the impurity solver and then determine Ud, by hnding the value that 
best matches the full spin-polarized MQSGW spectra. The procedure starts from the non- 
spin-polarized MQSGW band structure without magnetic long-range order. We then allow 
spontaneous magnetic long-range order by embedding a polarized impurity self-energy for 
the Cn-3d electrons computed in a local GW approximation. We hnd that indeed magnetic 
ordering associated with Cu-3d is captured by spin-polarized local MQSGW using a static 
value of Ud and Jh and spectral properties such as energy gap are very similar in value 
to the full spin-polarized MQSGW spectra. In Fig. 2(b), we allowed Ud to vary between 
8-13 eV (at hxed Jh = 1.4 eV) and we plot the size of the indirect gap. The gap size of this 
method matches the gap of spin-polarized MQSGW when Ud ~ 12 eV. If this choice of Ud 
and Jh is correct, the resulting spectra must be similar to the prediction of spin-polarized 
MQSGW method. We show this comparison in Fig. 2(c) to conhrm a good match. In 
addition, the relative position of the Gu-d band (the lowest energy conduction band at S) 
to the La-d band (the lowest energy conduction band at Y) is also well matched justifying 
the approximation of S^‘"(ia;„) ~ (icun, = 0). T,^^{iujn = 0) for Cu-dx 2 _y 2 orbital differs 
from nominal double counting energy 3^ by only 1%, highlighting again the advantages of 
using a broad window and narrow orbitals. 

We now discuss the magnetic moment associated with Gu and the electronic excita¬ 
tion spectra of La 2 Gu 04 by using MQSGW-I-DMFT (with Ud = 12.0eR, Jh = 1.4eV) in 
which the impurity is solved by the numerically exact GTQMG 3, 171 and compare them 
with other methods. LSDA does not have a magnetic solution. In contrast, spin-polarized 
MQSGW, QSGW 0, and MQSGW+DMFT predict 0.7 ps, 0.7 /ub, and 0.8 /ub, respec¬ 


tively. This is consistent with 


large range O.d/xs ~ 0.8/1^ 37-39] 


experimental measurements, although the later span quite 


In the low-energy spectrum of La 2 Gu 04 , LSDA does not have an insulating solution; there 
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is a single non-magnetic solution with zero energy gap as shown in the bandstructure(Fig. 
3(a)) and total density of states (Fig. 4(a)). The non-spin-polarized MQSGW also predicts 
metal as shown in Fig. 4(a), but the two bands of primarily Cu-d^^-y^ character near the 
Fermi level are well-separated from the rest of the bands (dashed lines in Fig. 3(b)). Spin- 
polarized MQSGW calculation (dashed lines in Fig. 3(c)) yields qualitative different results 
from LSDA and non-spin-polarized MQSGW calculation. The two Cu-dx 2 -y 2 bands are now 


well separated from each other with a bandgap of 3.4 eV. Spin-polarized QSGW [27| also 
yields insulating phase with a gap of 4.0 eV. In the experiment, the larger direct gap, as 


measured by optics, is ~ 2eV 


40, 




We show that these dehciencies of LDA, QSGW and MQSGW in the low-energy spectra 
can be remedied by adding all local Feynman diagrams for the Cn-d orbitals using the DMFT. 
The LDA-I-DMFT calculation in Fig. 4(a), carried out by the all-electron LDA-I-DMFT 


method 


31 


35|, predicts reasonable gap of 1.5 eV and 1.8 eV in PM and AFM phases, in 


good agreement with experiment and previous LDA-I-DMFT studies 


31 


42 


451. Within 


MQSGW-I-DMFT, we find gaps of 1.5 eV and 1.6 eV in PM and AFM phases, respec¬ 
tively, as shown in Fig. 4(b). The excitation spectra of MQSGW-I-DMFT in PM and AFM 
phase as shown in Fig. 3(b) and 3(c) are very similar as both are insulating with well 
separated Cu-dx 2 -y 2 bands, which is now also substantially broadened due to large scatter¬ 
ing rate in Hubbard-like bands. In addition, MQSGW-I-DMFT improves the line-shape of 
LDA-I-DMFT. Near the top of the valence bands with oxygen p character, the lineshape 
within LDA-I-DMFT is too sharp in comparison to the experiments as shown in Fig. 4(c). 
By treating oxygen p levels within GW, the lineshape becomes smoother and in a better 
agreement with experiments. 

In the high energy region of La 2 Gu 04 , the most distinctive difference is the position of 
La-/ peak. It appears at ~ 3eV within LDA and LDA-I-DMFT, but at around ~ 9eV, in the 
inverse-photoemission spectra (cyan dotted line in Fig. 4(a)) 36j. By treating La-/ within 
GW approximation, it appears at ~ lOeV within MQSGW and MQSGW-I-DMFT. The 
underestimation of unoccupied La-/ excitation energy is attributed to the local approxima¬ 
tion to the electron self-energy within LDA. Within LDA, Hartree and exchange-correlation 
potential applied to La-/ orbitals are orbital-independent since charge density is averaged 
over 14 different m channels In contrast, these potentials within MQSGW are orbital- 


dependent and non-local. The effect of orbital-dependent potential can be tested within 
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LDA+U approaches, since LDA+U adds orbital-dependent potential and subtracts orbital- 
independent potential explicitly j^. From LDA-I-U approaches, we can also understand 
MQSGW better since LDA-I-U can be regarded as a local and static approximation to GW 
approximation j^. According to M.T.Gzyzyk and G.A.Sawatzky ^|, La-/ peaks shift from 
Ei7’-|-3eV to Ei?-|-3eV-|-U/2 with U=lleV for La-/. 

We also tested our proposed scheme with one more charge transfer insulator, NiO. Fig. 
5(a) shows the frequency dependence of Ud and Wd for the Ni-dd orbitals in the low-energy 
region. In contrast to Wd, Ud is almost constant up to 5eV. Static Ud is 9.6 eV. In the 
high energy limit, Ud and Wd approach the bare value of 26.0 eV. Galculated Jh for the 
Ni-3d orbitals has negligible frequency dependence and static Jh is 1.4eV. Fig. 5(b) shows 
the total density of states of NiO within LDA-I-DMFT and MQSGW-I-DMFT in its para¬ 
magnetic phase. Photoemission/inverse photoemission data are also plotted for comparison 


^|. The LDA-I-DMFT calculation is being carried out by the all-electron LDA-I-DMFT 


method 


35j with Ud = lOeV, Jh = 0.9eU and nominal double-counting energy. In the 


paramagnetic phase, LDA-I-DMFT and MQ SGW-I-DMFT predict insulator in an agreement 
with previous LDA-I-DMFT studies 49, 50|, but MQSGW-I-DMFT improves the line-shape 
of LDA-I-DMFT. Near the top of the valence bands, the lineshape within LDA-I-DMFT 
is too sharp in comparison to the experiments. By treating oxygen p levels within GW, 
the lineshape becomes smoother and in a better agreement with experiments. In the an¬ 


tiferromagnetic phase, magnetic moment associated with Ni-d orbital 
MQSGW-I-DMFT, in agreement with experimental value of 1.6-1.9 ps 


s is 1.6 aB within 


49|,y,l52|. 


In summary, we introduced a new methodology within MQSGW-I-DMFT and tested 
it in the classic charge transfer insulator La 2 Gu 04 and NiO. Our methodology predicts a 
Mott-insulating gap in the PM phase, thus overcoming the limitation of LDA and QSGW. It 
yields more precise peak positions of the La-/ states in La 2 Gu 04 and valence band lineshape, 
thus improving the results of LDA-I-DMFT. The method should be useful in understanding 
electronic excitation spectrum of other strongly-correlated materials, in particular, those 
where precise position of both the itinerant and correlated states is important. 

Methods 

Our approach is carried it out entirely on the Matsubara axis, which requires a different 
approach to the quasiparticle self-consistency in GW |53|, called Matsubara Quasiparticle 
Self-consistent GW (MQSGW), where the quasiparticle Hamiltonian is constructed by lin- 
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earizing the self-energy and renormalization factor 5^. Working on the Matsubara axis 
is numerically very stable, provide a natural interface with advanced DMFT solvers such 
as continuous-time quantum Monte-Carlo (CTQMC) [^, 17| and has very good scaling in 
system size as in the space-time method, (see supplementary note on Matsubara QSGW 
calculations). 

For DMFT, it is essential to obtain bandstructures in a hne enough crystal momentum 
(k) mesh to attain desired frequency resolution of physical quantities. To achieve such 
momentum resolution, we use a Wannier-interpolated MQSGW bandstructure in a large 


energy window using Maximally localized Wannier function (MLWF) and 


structed local projector in a hne momentum mesh. In contrast to SrVOs 12l-ll5| where a 


hen con- 


set of t 2 g states is reasonably well separated from the other bands, correlated 3d orbitals in 
La 2 Cu 04 and NiO and are strongly hybridized with other itinerant bands. In this case, it 
is necessary to construct local projectors from states in a wide enough energy windows to 
make projectors localized near the correlated atoms. We constructed local projectors in the 
energy window Ep ± lOeV in which there are ~82 bands at each k point, where Ep is the 
Fermi level for La 2 Gu 04 . For NiO, we constructed local projectors in the energy window 
of Ep — lleV to Ep + lOeV. Then we conhrmed that absolute value of its overlap to the 
muffin-tin orbital (of which radial function is determined to maximize electron occupation 
in it) is more than 95%. Our choice of energy window is justihed by the Cu-3d spectra 
being entirely contained in this window. Using constructed MLWFs, we dehned our local- 
projector Pi,n(k) = where lURi(r) is MLWF with an index i, 

'0rik(r) is quasiparticle wavefunction with an index n, and W is the number of k points in 
the hrst Brillouin zone. 


Static Ud and Jp are evaluated by a modihcation of the constrained RPA method [56|, 


which avoids screening by the strongly hybridized bands. This screening by hybridization is 
included in our large energy window DMFT. For details, see supplementary note on Ud and 
Jp- We divide dynamic polarizability within MQSGW approximation xqp into two parts, 
Xqp = Xqp + Xq^- Here, Xqp is dehned by all transitions between the states in the energy 
window accounted for by the DMFT method {Ep ± lOeU for La 2 Gu 04 and Ep — lleV 
to Ep -|- lOeU for NiO). Using Xqi^j we evaluate partially-screened Goulomb interaction 


U ^{r,r',k,iun) = 
Slater’s integrals [57 


/ ^(r,r',k) — Xgp^(r, r', k, ia;„) and parametrize static Ud and Jp by 
58| . where V is bare Goulomb interaction. 
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The Feynman graphs included in both MQSGW and DMFT (double-counting) are the 
local Hartree and the local GW diagram. They are computed using the local projection 
of the MQSGW Green’s function (Gqp) GQp{iuJn) = '^]iPO^)GQp(k,i^n)PK^) ^^e 
local Goulomb matrix constructed from Slater’s integrals. For the details, see supplementary 
note on double counting energy. 
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Figure 1. Atomic structure and first Brillouin Zone of La 2 Cu 04 . (a) Atomic structure of La 2 Cu 04 
in the single face-centered orthorhombic phase. Lanthanum atoms are represented by green spheres, 
copper atoms by blue spheres in the blue octahedrons, and oxygen atoms by red spheres. The 
structure is characterized by an alternating rotation of successive CuOe octahedra along the x 
direction, (b) First Brillouin zone of single face-centered orthorhombic phase. Red lines show the 
path along which electronic bandstructures are plotted in Fig. 2(c) and Fig. 3. 
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Figure 2. Hubbard U associated with Cu-3d orbitals in La 2 Cu 04 . (a) Frequency dependence of 
Wd (dashed lines) and Ud (full lines) of La 2 Cu 04 with a Xq'p defined in the energy window Ep 
lOeH. Real and imaginary parts of the parameter are marked by red and blue colors, respectively, 
(b) Bandgap dependence on Ud, in La 2 Cu 04 , evaluated with impurity self-energy within spin- 
polarized GW approximation with J//=1.4eV. The Black dashed line represents bandgap within 
spin-polarized MQSGW. (c) Spectral function of La 2 Gu 04 with Ud=12eV and J//=1.4eV. The 
black dashed-lines show bandstructures within spin-polarized MQSGW 
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Figure 3. The low-energy spectral function of La 2 Cu 04 . (a) Electronic bandstructures of La 2 Cu 04 
within LSDA and spectral functions from (b) non-spin-polarized MQSGW-I-DMFT (c) and spin- 
polarized MQSGW-I-DMFT calculations along the path shown in Fig. 1(b). The Dashed lines in (b) 
and (c) represent electronic bandstructures within non-spin-polarized MQSGW and spin-polarized 
MQSGW, respectively 
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Figure 4. The density of states of La 2 Cu 04 . (a) Total density of states of La 2 Cu 04 from LDA 
(magenta), LDA+DMFT(green), MQSGW (red), and MQSGW+DMFT (blue). Full lines and 
dashed-lines represent quantities within non-spin-polarized and spin-polarized versions of each 
calculation, respectively. The cyan dotted line shows photoemission/inverse-photoemission data 


36l |. The Positions of La-/ peaks are marked by arrows, (b) A zoom-in view in the low-energy 


region, (c) The overlap of total density of s^a^es of La 2 Cu 04 within LDA-I-DMFT as well as 
MQSGW-I-DMFT and photoemission/inverse-photoemission data 
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Figure 5. Hubbard U associated with Ni-3d orbitals in NiO (a) Frequency dependence of Wd 
(dashed lines) and Ud (full lines) of NiO, with a Xqp defined in the energy window in Ep — lleV to 
Ep + lOeF. Real and imaginary parts of the parameter are marked by red and blue colors, respec¬ 
tively. (b) Total density of states of NiO within LDA+DMFT(green) and MQSGW-l-DMFT(blue). 


The cyan dotted line shows photoemission/inverse-photoemission data 481 


18 

















































arXiv:1504.07569v2 [cond-mat.str-el] 29 Aug 2016 


Supplementary information for First-principles treatment of 
Mott insulators: linearized QSGW+DMFT approach 

Sangkook Choi/’^ Andrey Kutepov,^ Kristjan Haule/’^ 

Mark van Schilfgaarde,^ and Gabriel Kotliar^’^ 

^ DMFT-MatDeLab center, Upton, New York 11973, USA 
“^Department of Physics and Astronomy, 

Rutgers University, Piscataway, New Jersey 08854, USA 
^Department of Physics, Kings College London, 

Strand, London WC2R 2LS, United Kingdom 


1 


Computational details: basis set 


Calculations are performed using our relativistic ^in-polarized, full-potential, linearized- 
augmented-plane-wave package (RSPFLAPW) |l|, y], which is based on full-potential lin¬ 
earized augmented plane wave plus local orbital method. For La 2 Cu 04 , experimental lattice 
constants and atomic positions at the low-temperature orthorhombic phase are used. The 
following parameters for the basis are used: muffin-tin (MT) radii (Rmt) in Bohr radius are 
2.2 for Lanthanum(La), 2.1 for Cu, and 1.4 for Oxygen(O). Wave functions are expanded by 
spherical harmonics with / up to 4 for La, 4 for Cu, and 3 for O in the MT spheres, and by 
plane waves with the energy cutoff determined by Rmt,Cu x Kmax = 7.8 in the interstitial 
(IS) region. The Brillioun zone was sampled with 3x3x3 k-point grid. For the density 
functional theory calculation, local spin density approximation (LSDA) is employed. For the 
GW calculation. The following convergence parameters are used. Product basis is expanded 
by spherical harmonics with I up to Imax =5 in the MT spheres and Rmt,Cu x K^ax = 9.0 in 
IS region. Unoccupied states with an energy up to 500 eV from Fermi energy are taken into 
account for polarizability and self-energy calculation. The projector is constructed using 
MLWF and interpolated in a 10 x 10 x 10 k-grid. Spin-orbit coupling was not included. 



For NiO in the rocksalt structure with a lattice constant 


following 


parameters for the basis are used: muffin-tin (MT) radii {Rmt) in Bohr radius are 2.12 


for Nickel(Ni) and 1.77 for O. Wave functions are expanded by spherical harmonics with I 
up to 6 for Ni and 6 for O in the MT spheres, and by plane waves with the energy cutoff 


determined by Rmt^ni x Kmax = 6.7 in the interstitial (IS) region The Brillioun zone was 
sampled with 5x5x5 k-point grid. For the GW calculation. The following convergence 
parameters are used. Product basis are expanded by spherical harmonics with I up to l^ax 
=6 in the MT spheres and Rmt,Cu x Kmax = 10.0 in IS region. Unoccupied states with 
an energy up to 680 eV from Fermi energy are taken into account for polarizability and 
self-energy calculation. The projector is constructed using MLWF and interpolated in a 
15 X 15 X 15 k-grid. Spin-orbit coupling was not included. 

Matsubara QSGW calculations 

The electron self-energy can be systematically expanded in terms of the dressed Green’s 
function G and the screened Coulomb interaction W. Within GW approximation, we keep 
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only the first term of the series expansion of self-energy in W: 

Sciy(r, r', k, s, iun) = - ^ f dTG{r, r', R, s, t)W{t, r', R, (SI) 

R d 

where r is the position vector in a unit cell, k is the crystal momentum, R is the lattice vec¬ 
tor, Un is Matsubara frequency, r is Matsubara time, and s is spin index. Within Matsubara 
quasiparticle self-consistent GW (MQSGW) approximation, we calculate the dynamic polar¬ 
izability and the electron self-energy by using quasiparticle Green’s function {Gqp) instead 
of Full GW Green’s function. First, we construct quasiparticle (QP) Green’s function using 
the quasiparticle Hamiltonian {Hqp) by 


Ggpir, r', k, iujn) = iujJir - r') - HQp{r, r', k). 


(S2) 


For the hrst iteration, we regard Hamiltonian within local density approximation as Hqp. 
Next, the screened Goulomb interaction Wqp is evaluated using the dynamic polarizabil¬ 
ity XQp{r,r',k,' iuJn) by WQp{r,T' ,k,iuJn) = 'F"^(r,r',k) - XQp{r,r' ,k,iuJn) within random 
phase approximation (RPA): 


XQp{r,r',k,'iuJn) = EE/ drGgp/r, r', R, s, r/Ggp/r', r, -R, s, -r) 




s R 


where V is bare Goulomb interaction. Then, we calculate MQSGW electron self-energy 
(Egp) by 

EQp(r, r', k, s, iun) = - ^ j drGQpir, r', R, s, t)Wqp{y, r', R, 

R d 

Then, we constructed quasiparticle Hamiltonian with linearized self-energy and renormal¬ 
ization factor, E“^(r,r',k) = 1 — SSgp/r, r', k, = Q)/d{iuJn) by 


Hqpik) = Z^/\k) (Rp(k) + EQp(k,za;„ = 0)) Z^/\k), 


(S5) 


where Rp(k) is Hartree Hamiltonian. This process is repeated until self-consistency is 
attained. 


Ud and Jh 

The static Ud and Jp are evaluated by a method similar to constrained RPA , but here 
we avoid screening by the strongly hybridized bands, which is included in our large energy 
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window DMFT. We divide dynamic polarizability xqp iiito two parts, xqp = Xqp + Xqp^- 
Here, Xq'p is defined by all transitions between the states in the energy window acconnted 
for by the DMFT method (Ep ± lOeH for La 2 Cn 04 and Ep — lleV to Ep + lOeV for NiO): 


unocc occ 

in the window in the window 


xSp(r,r',k,ia;„) =-2^ 


k' 


V'nk'(r)V'mk'+k(r)V'nk'(r')V'mk'+k(r')- 


~ (En-k' — -Enk'+k)^ 


(S 6 ) 


Using Xqp^, we evalnate partially-screened Coulomb interaction by U~^{r,r','k,iun) = 
U“^(r, r',k) — XQp^(r^', k, icun)- Then, we calculate U = Eq and Jp = {E 2 + F4)/14 by 
using Slater integrals [Gj, 


F, = 


dvr 


Ck2k + 1 


^ ^ (h 2 mi I (h 2 m, 2 ^^m 2 —m 3 | ^2mz) 


mi,m 2 ,m 3 ,7714 


(ST) 




where HR, 4 (r) is Wannier function centered at R with index m and Yim is spherical har¬ 
monics. Here, Cq = 25, C 2 = 20/49 and C 4 = 100/441. 


DMFT: local projection and embedding 


With non-spin-polarized MQSGW Hamiltonian (Hqp), we solve DMFT self-consistent 
equation. 


^ ^ Ei,ni (k) HQp(\i, iUJji^ Semf)e(i(k, S, ZCUfi) 

k 

Eimp Aj,7jp(s, ibj'^ Sjmp(s, 


-1 


PlnS^) 


711,712 


(S 8 ) 


Here, Pi,n(k) = Yhp < hFaj|'^„k > e®^’^/\/^ is projector to the correlated subspace at each 
k and 'i/’nk(r) is quasiparticle wavefunction with an index n. is the number of k points in 
the first Brillouin zone. Eimp and Ajmp are impurity level energy and hybridization function, 
which are inputs to impurity solver, f^embedded = Fl(k) (f:imp{s,iujn) - F(k) 

is embedded self-energy with impurity self-energy (Sjmp) and double-counting correction 
(Sdc)- 
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Double counting energy 


The Feynman graphs included in both MQSGW and DMFT (double-counting) are the 
local Hartree and the local GW diagram. They are computed using the local projection of 
the MQSGW Green’s function (Gqp) &Qp{iuJn) = -^'^y_P(k)GQp(k,icdn)PH^) and the 
local Goulomb matrix Uikij'- 

= 0-)Uia, - f (S9) 

k,l=Cu-d k,l=Cu-d 

where, Wl^^i{iUr,)=Uikji + Y.n^npq=Cu-dUininiP^pqni^u}r^)W^^^^^^^ and P^pgn{i<^n) = ‘^ J dv 
^QP,n,pi'^)^QP,q,m{~pG‘^^. The Goulomb matrix Uikij is constructed by using Slater in¬ 
tegrals of Fq = 12.0eG, F 2 = 12.1eG, and F 4 = 7.5eV for La 2 Gu 04 and Fq = 9.6eG, 
F 2 = 12.1eG, and F 4 = 7.5eV for NiO in the following way: 


F' 


ml,m2,m3,m4 , 


iPn) = 5'mi m' *5^, m' *5' , S 

iLj m2,m2 7713,7713 7774, 
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E 

A:=0 


dvr 


2k+ 1 


{y2m[\ ) (^2m2^fcgl y2mp F {iVn) 


(SIO) 


where S is the transformation matrix from spherical harmonics to cubic harmonics. Finally, 
for the stable numerics, we approximated = 0 ) since these low order 

diagrams are dominated by the Hartree-Fock contribution. 
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